A dynamical symmetry of the spherical dust collapse 
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By linearly scaling the initial data set (mass and kinetic energy functions) together with the initial 
area radius of a collapsing dust sphere, we find a symmetry of the collapse dynamics. That is, the 
linear transformation defines an equivalence class of data sets which lead to the same end result as 
well as its evolution all through. In particular, the density and shear remain invariant initially as well 
as during the collapse. What the transformation is exhibiting is an interesting scaling relationship 
between mass, kinetic energy and the size of the collapsing dust sphere. 

PACS numbers: 04.20 Dw, 04.20 Ex 



In classical general relativity, one of the key outstanding questions is the end result of a 
gravitationally collapsing object, does it end in a black hole (BH) or in a naked singularity 
(NS) [1, 2] ? There exists no method to address this question in general terms, and hence 
one is left to study various examples of collapsing systems [3] with a view to gain some 
insight. Recently, the most pertinent question, what is the physical process that causes NS, 
has also been addressed for spherical dust as well as for the general Type I matter field 
[4,5]. The answer is that it is the shearing force in the collapsing cloud which is responsible 
for NS. It has been proven that shear is necessary for NS while it has to be strong enough 
(beyond certain threshold value) to lead to NS. 

Recently, again for a spherical dust collapse, a linear transformation has been found 
which keeps density and shear of the collapsing sphere invariant for all through the collapse 
but for the initial instant [6]. That is the initial density and shear are not invariant under 
this transformation. In this letter we wish to complete the invariance for initial event as 
well by bringing in the radius of the sphere as well to scale linearly. It is now a symmetry 
of the collapse dynamics and we have an equivalent set of initial data leading to the same 
end result as well as its evolution all through the collapse. 

The inhomogeneous spherically symmetric dust collapse in n dimensions can be described 
by the higher dimensional Tolman-Bondi-Lemaitre (TBL) line element given by 

ds 2 = -dt 2 + - —dr 2 + R 2 dQ 2 n _ 2 , (1) 

1 + f(r) " 2 w 

where the area radius, R — R(t, r) and f(r) > —1. The evolution equation for R is 

R 2 -^ + m, (2) 

where the arbitrary function F(r) > is related to the matter density as 



. , (n-2)F'(r) 



The measure of anisotropy is given by the shear scalar 
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n-2 fit' R\ I n-2 [{RF 1 - (n - 1)R'F} + R n - 3 (Rf - 2R'f)] 



V2(n-1)\# R) ys(n-l) R^ R' (F + fR^) 1/2 

(4) 

Now using the initial condition R = r at t = t% (initial epoch) the expressions for initial 
density and shear are 

Pl (r) = p(U,r) = ^r 2 - n F'(r) (5) 

and 



n-2 [{tF' -(n-l)F} + r"- 3 (rf -2/)] 

ai = <j{ti,r) = * — — — — - (6) 

y8(n-l) (F + /r"- 3 ) 1/2 

For the finiteness of the initial density and shear at the centre, the regular functions F 
and / must have the following series expansion 

F{r) = For 11 - 1 + Fir" + F 2 r n+1 + (F ^ 0) 

(7) 

f(r) = for 2 + fir 3 + f 2 r* + (/„ ^ 0) 

For smooth initial density profile we write 

Pi(r) = Po +Pir + p 2 r 2 + (8) 

So using these series expansions in equation (5) we have the relation among the different 
coefficients as 

_ ( n +j-l)(n-2) 
Pj — 7, b 3i 1,2,... 



The evolution equation (2) can be integrated to give 



t-U = 



{n-l)y/F 



1 fr n - 3 , n=i 1 fR n - 3 



r * 2 F 1 [-,a,a+l, — ] - R * 2 F 1 [-, a, a + 1, - 



2' ' ' F *-*l 2 >~>- ■ -> F 



(9) 



where 2 Fi is the usual hypergeometric function with a = h + — Kr. 



If t = t s (r) denotes the time of collapse of the r-th shell to the singularity then 
R(t s (r),r) = and we have from the above equation 

2 n-i 1 f r n-3 

ts(r) ~ U = (n-l)VF r ^ 2Fl[ 2 ' °' ° + (10) 

The formation of trapped surface is signaled by R n ~ 3 = F(r) at t = t ah (r) and so we 
write, 

2 r i T i 1 fr n ~ 3 2F^ 1 

t ah {r)-U = - --7=2^1 -,a,o+l,-i-=- 2 Fr -,a,a+l,-/ (11) 

(n-l)VF 2 i< rc-1 2 

Using the series expansion (7) if we proceed to the limit as r — > in equation (10) then 
we have the limit of central singularity as 
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to — ti + 



(n-iyrs 



2^1^,0,0+1,-^] 



(12) 



Hence the time difference between the formation of apparent horizon and central singu- 
larity is given by (upto leading order in r) 



tah - to 



' f n — 3 



{n-l)y/F 



— 2 F 1 [-,a,a+l,-—] 



(n-1) f (h FA 3 ^ 9 / . 



This reduces in 4 and 5 dimensions respectively to 



(13) 



tah — to — — 




2(/ F 1 - hFp) 3 5 7 _^ 
+ ^5/2 2 ^ U 2'2'2' i^o 



0(r 2 ) 



(14) 



tah — to 



(2/1F0-/0F1) , (/ Fi - / /i - 2/iJb) 



+ 



0(r 2 ) 



(15) 



Obviously it is NS when t a h > h and BH otherwise. It is the data set {F, /} that 
determines the ultimate end result in terms of BH or NS for the collapsing sphere. 

Mena etal [6] have proposed the following linear transformation on the data set I = {F, /}, 



{F,f}^{a^F,af} 



(16) 



(a > a constant) under which both p,a remain invariant but not initially, i.e pi, at are 
not invariant. They have however noted that it is this discontinuous behaviour which 
is responsible for having NS for initially vanishing shear or density contrast. It is this 
discontinuity which compensates for absence of initial density contrast or shear. It is the 
reflection of the fact that it is always possible to generate an NS for any arbitrary density 
profile by suitably choosing the energy function / [7]. It is however interesting that this 
transformation also picks up this general aspect. 

Our concern here is to seek a symmetry of the collapse dynamics - a true invariance. It 
stands to reason that when we linearly scale mass and energy functions, we should also scale 
the radius R of the sphere to keep everything in balance. The above transformation should 
thus be supplemented by 



R-> yfaR. 



(17) 



Now this is the symmetry of the complete collapse dynamics. It can be easily verified 
that the density and shear both initially as well as otherwise do remain unaltered under 
the transformation. Also the time difference (t a h — to) remains invariant indicating the 
complete symmetry of the collapse process. For marginally bound collapse we can set / = 



4 



in the transformation without disturbing the invariance. 

We have thus completed the linear transformation proposed by Mena etal [6] to make 
it as a symmetry of the collapse dynamics. What it indicates is that if we scale the 
radius R by a, then mass function F will naturally scale as a n ~ l and the energy / as a 2 
(dimensionally F — > RS n ~^F). The transformation is therefore physically quite transparent 
but the remarkable fact is that the overall dynamics also respects it all through. This is 
the dynamical symmetry of spherical collapse. It is though established for the dust collapse 
but it should be true in general for other matter fields as well. 
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